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Abstract We review the basic physics behind light in- 
teraction with plasmonic nanoparticles. The theoretical 
foundations of light scattering on one metallic parti- 
cle (a plasmonic monomer) and two interacting parti- 
cles (a plasmonic dimer) are systematically investigated. 
Expressions for effective particle susceptibility (polar- 
izability) are derived, and applications of these results 
to plasmonic nanoantennas are outlined. In the long- 
wavelength limit, the effective macroscopic parameters 
of an array of plasmonic dimers are calculated. These 
parameters are attributable to an effective medium cor- 
responding to a dilute arrangement of nanoparticles, i.e., 
a metamaterial where plasmonic monomers or dimers 
have the function of "meta-atoms". It is shown that pla- 
nar dimers consisting of rod-like particles generally pos- 
sess elliptical dichroism and function as atoms for planar 
chiral metamaterials. The fabricational simplicity of the 
proposed rod-dimer geometry can be used in the design 
of more cost-effective chiral metamaterials in the optical 
domain. 
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1 Introduction 

Light waves cause the charged constituents of matter, 
electrons and nuclei, to oscillate. These moving charges 
in turn emit secondary light waves, which interfere with 
the incident wave and with each other. By an appro- 
priate spatial averaging over the fields at atomic length 
scales, the collective response of all constituents can be 
derived. This macroscopic response strongly depends on 
the type of material. This is how the conventional mate- 
rial equations that relate the strength and displacement 
fields in the electromagnetic wave are introduced, and 
this is how electrodynamics and optics of homogeneous 
media is formulated. 

Along similar lines of reasoning, one can engineer the 
building blocks ("meta-atoms") with feature sizes smaller 
than the wavelength of light, and still perform a spatial 
averaging over the fields at the "meta-atom" level. This 
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Figure 1 (Color online) (a-b) Schematic illustration and (c- 
d) plasmonic dimer representation of 3D and 2D enantiomeric 
meta-atoms. 



results in "macroscopic" effective parameters for such ar- 
tificial composite materials. The ability to design the 
meta-atoms in a largely arbitrary fashion adds a new 
degree of freedom in material engineering, allowing to 
create artificial materials (metamaterials) with unusual 
physical phenomena rare or absent in nature. Examples 
include media with negative refractive index [l|, Q , hy- 
perbolic or indefinite media Q , and inhomogeneous me- 
dia capable of guiding the radiation around objects ("op- 
tical cloaking" @, S@, 0|) or m t° objects ("optical black 
holes" Q) using the concept of transformation optics 

Another group of unusual physical phenomena in meta- 
materials results from their ability to transform polar- 
ization of light. One example is giant optical activity (Io| 
in composite materials containing spiral-like or otherwise 
twisted meta-atoms (see, e.g., figure QJi). Recently, pla- 
nar chiral metamaterials (PCMs) were introduced [111 . 
lit 13 1 where meta-atoms possess two-dimensional (2D) 
rather than three-dimensional (3D) enantiomeric asym- 
metry (figure [Tp) . While conventional chiral PCMs sim- 
ply resemble naturally occurring gyrotopic media (e.g., 
optically active liquids), PCMs are distinct from both 
3D chiral and Faraday media in that their polarization 
eigenstates are co-rotating elliptical rather than counter- 
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rotating elliptical or circular [l2, 1^, 14 1. This leads to ex- 
otic polarization properties, e.g., asymmetry in transmis- 
sion for a forward- vs. backward-propagating circularly 
polarized incident wave, without nonreciprocity present 
in Faraday media. Such a rich variety of polarization 
properties in compact-sized planar structures make chi- 
ral metamaterials promising for polarization sensitive in- 
tegrated optics applications. 

Taken generally, the analysis of artificial materials 
is a complex theoretical and computational problem. In 
the theoretical research of metamaterials, two areas can 
roughly be identified. On the one hand, there is the ac- 
tive research concerning various approaches to metama- 
terial homogenization, i.e., the derivation of macroscopic 
material parameters (in particular, permittivity and per- 
meability tensors) out of the electromagnetic properties 
of a single meta-atom 15, 0,03, El]- On the other hand, 
a lot of effort is invested in the study of light interac- 
tion with single nanoparticles (see [19j and references 
therein) . Metallic (plasmonic) nanoparticles attract spe- 
cial attention [20| because they can exhibit strong reso- 
nances for the wavelengths much larger than their fea- 
ture sizes. 

Both experimentally and theoretically, one typically 
takes advantage of the scalability of Maxwell's equa- 
tions to consider structures with millimeter-scale feature 
sizes in the microwave and radio-frequency (RF) regime 
(HI . [l2l . [T^ . Not only is the metamaterial fabrication 
on these length scales much easier, but also their theo- 
retical analysis - meta-atoms have a close resemblance 
to antennas in their ability to connect between propa- 
gating radiation and localized near fields. The extensive 
knowledge in the field of antenna engineering can then 
be tapped into. Semi-analytical models exist [2lj that 
offer insight into field interaction processes in the an- 
tennas, and numerical tools for practical antenna design 
and optimization are well established (22| . 

The optical wavelength range, which is where func- 
tional metamaterials attract the most interest, remains 
far more challenging. At optical frequencies, the metal 
can no longer be treated as a perfect electric conductor. 
Because of this fundmental difference the semi-analytical 
models used in the microwave regime can not be used 
in the visible and near-infrared spectral range, and the 
demands in computational resources increase tremen- 
dously. Fabrication of metamaterials in the optical do- 
main also remains a challenge because of the small fea- 
ture sizes (about 10 to 100 nm). 

These difficulties notwithstanding, recent progress in 
nanotechnology have enabled the fabrication of optical 
antennas (nanoantennas) (23l 24 1 and opened many ex- 



citing possibilities towards nanoantenna applications. For 
example, it has been recently demonstrated that nanoan- 
tennas can enhance [25[ and direct [26j the emission of 
single molecules and that they can play a key role in 
sensing applications [27j. Great potential in improving 
the efficiency of solar cells should also be mentioned [28]. 



Such a combination of promises and challenges in op- 
tical metamaterials creates a strong need for a system- 
atic, ab ibitio, theoretical description of metamaterials, 
which, on the one hand, would be general enough to 
span the range of applicable frequencies from RF to op- 
tical, and on the other hand, would be simple enough to 
facilitate the direct prediction of metamaterial proper- 
ties based on their geometry without having to resort to 
full-scale numerical simulation in time-consuming adap- 
tive optimization techniques. Such a theoretical descrip- 
tion would also be beneficial in order to be able to sim- 
plify the geometry of functional (e.g. chiral) meta-atoms 
so as to facilitate their cost-effective fabrication in the 
telecommunication and optical domains. 

In this paper, we formulate such a theoretical de- 
scription by investigating the theoretical foundations of 
light interaction with plasmonic nanoparticles. We start 
with one metallic particle (a plasmonic monomer) in an 
external electromagnetic field, and consider the general 
expression for its susceptibility in terms of the formal 
Green's function-based solution of the Maxwell equa- 
tions. We then set out to simplify this general solution 
for a nanoparticle of a simple geometry in the approx- 
imation of its size being much smaller than the wave- 
length (the Rayleigh approximation), moving on to the 
corrections to this approach accounting for the retarda- 
tion effects. Further we introduce a new semianalytical 
method to solve the problem of electromagnetic wave 
scattering on metallic cylindrical wire, which is both ac- 
curate and computationally efficient. The results turn 
out to be equally applicable in the microwave and opti- 
cal frequency range and can be used for improved designs 
for plasmonic nanoantennas. 

We then apply a similar approach to a pair of closely 
located nanoparticles (a plasmonic dimer). By account- 
ing for the interaction between the partices in the dimer, 
we arrive at analytic expressions of the dimer's polar- 
izability. It is then possible to determine the effective 
macroscopic parameters of a dilute medium that con- 
tains the corresponding dimers as "meta-atoms". It is 
shown that dimers of planar geometry without an in- 
plane mirror symmetry (figure[T}i) exhibit crystallographic 
and spectral properties of planar chiral metamaterials. 
It is concluded that a geometry as simple as an arrange- 
ment of nanorods can support functional metamaterial 
properties that were previously reported in more intri- 
cately shaped meta-atoms. 

The structure of the paper is as follows. In Section 
[2j the systematic approach to single plasmonic parti- 
cles (monomers) is presented. We start from the general 
Green's function solution and review the computational 
strategies for numerical and semianalytical approaches, 
ending up with the expressions of effective particle sus- 
ceptibility (polarizability) . In Section [3J a similar ap- 
proach is applied to plasmonic dimers. The effective po- 
larizability of a dimer is determined and the expressions 
for effective macroscopic properties such as effective per- 
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mittivity and permeability is derived. Application of the 
developed formalism to 2D planar chiral double-rod plas- 
monic dimers is then presented. Finally, Section [4] sum- 
marizes the paper. 



2 Plasmonic monomers 

2.1 General solution 

We describe a nanoparticle by a complex, frequency de- 
pendent relative dielectric permittivity e r . Electric field 
E (r) in the presence of such a nanoparticle is given by 
a solution of the homogeneous Helmholtz equation [29j 



V x V x E (r) - k 2 e (r) E (r) = 0. 



(1) 



Here = w/c is a free space wave number and e (r) is 
equal to e r inside the volume V occupied by the particle 
and to 1 otherwise. Representing the total field E (r) = 
E* (r) + E s (r) as a sum of incident E J (r) and scattered 
E s (r) fields and talking into account that the incident 
field has to satisfy a free-space Helmholtz equation 



V x V x E 4 (r) — fc 2 E l (r) = 



(2) 



one can reformulate the scattering problem ([T]) in the 
form of the inhomogeneous Helmholtz equation [29j 

VxVxE s (r)-fc 2 E s (r)=i^ f(r) (3) 

with an equivalent current density 

r(r) = -i W £ Z\e(r)E(r), (4) 

where Ae(r) = e (r) — 1. In this formulation one can 
clearly see, that the field scattered from a nanoparticle 
is generated by a current density induced in the particle 
by the total field E(r). 

Taking into account the explicit form of the induced 
current density ^ a general solution of the inhomoge- 
neous Helmholtz equation ^ can be written as a self- 
consistent integro-differential equation of the Lippman- 
Schwinger type [29( 



To find a general closed- form solution of the Lippman- 
Schwinger equation ([3]) we follow here the approach pre- 
sented in (30 . l3~0 | . It is well known that one can build 
a solution of equation ([5]) in an iterative manner by re- 
placing a self-consistent field in the integral by the field 
obtained on the left-hand side of the equation at the pre- 
vious iteration. Starting with the incident field one ends 
up with an infinite Born series expansion 32|. which in 
operator notations reads 



E = E l 



Y 



<rr>-o 



E 



'+( 



E' 



¥ 



E 4 



Y 2 



Ye 1 



Y-Y 



+7* 



tE' 



(8) 



In the last step we have used the formal summation of 
the operator power series. 

It is important to stress here that the Born expan- 
sion of a finite order can only be used in the case of rela- 
tively weak scattering. In the case of resonant scattering, 
which is of major importance in nanoplasmonics, finite- 
order Born approximation will lead to inaccurate results 
no matter how many terms in the expansion are taken 
into account (33J. At the same time, an infinite Born sc- 
ries provides an exact solution of the Lippman-Schwinger 
problem ([5]) both in the resonant and the non-resonant 
case. 

The main idea behind finding an exact solution of 
10 is to assume that the scattered field W (r) can be 
obtained via some dressed integro-differential operator 
acting on the incident field E' (r) directly [3CJ, [31| 



E*(r) = E(r! - E'(r) = lr ( ¥ + <g> V 



/ dV 5 (r,r')Y(r')/k(r')E*(r') (9) 
Jv 



or in operator notations 



E(r) =E 1 (r) + k z I 



<r> 1 



k 2 



V® V 



f dVj(r,r')4£(r')E(r'), (5) 
Jv 



where I denotes the three-dimensional unit tensor, (g> 
is the tensor product defined by (a ® b),^ = a,ibj and 



9(r,r') = 



Ak\r- 



4tt r — r' 



(6) 



E = E' 



TTE J 



(10) 



Here the unknown tensor a (r) has to be determined. 
Comparing relation (fTT)|) with the last line of relation 
JSJ) it can be immediately recognized that the dressed 
operator T must satisfy the Dyson-like equation 



YY. 



(ii) 



is the scalar Green's function. Denoting the integro-differential To find an unknown tensor a (r) we act with the 
^ ■ - " - operator (jTTJ) on an arbitrary plane harmonic incident 

field E 1 (r) = Eo exp (ik • r) , where k is the wave vector 
and Eo the amplitude of the incident field. This leads af- 
E = E J + t E. (7) ter a few simple steps to the following integro-differential 



operator in ([5]) by t the Lippman-Schwinger equation 
can be formally casted in operator notations as 
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equation for the unknown tensor a 



lr [ Y + ® V j y dV 



| ( Y (r') .9 (r, r') As (r') - g (r, r') (r')) c ikr '- 
Y(r')g(r,r')^e(r')fc 2 ( Y + ® V') | dV 



(r',r")Z\£(r")e ikr "f>E =0. (12) 



Since the amplitude of the incident field Eo can be 
arbitrary, the integral in this equation can be zero only 
for a vanishing integrand. That means that the expres- 
sion in the curly brackets should be equal to zero leading 
to 



Y(r')-Y-Y(r')fc 2 (Y + iv'sv 



d 3 r"g (r', r") As (r") e ik '( r "- r ') = 0. (13) 

v 

This equation can be easily solved to find the unknown 
tensor 

Y (r ) = (Y-Y(r))- 1 

where 



(14) 



Yw^fc 2 fY + iv®v) ^d 3 

3 (r,r')Z\e(r')e ik ( r '- r ) (15) 

plays the role of a self-energy operator. With tensor 
a (r) found, equation ([9]) gives an exact closed-form 
solution of the Lippman-Schwinger equation ([5]). 

It is instructive to reformulate this solution slightly. 
Recalling that Ae (r) = e (r) — 1 and e (r) = x(r) + 1, 
where x ( r ) is the electric susceptibility, one can rewrite 
the exact solution ([9]) in the following form 

E(r) = E l (r) + k 2 ^ Y + ® VJ J dV 

ff (r,r')Y(r')E*(r'), (16) 

where an effective susceptibility tensor of a nanoparticle 

X (r) = x ( r ) " ( r ) was introduced. Please note that 
such a general solution is a solution of the inhomoge- 
neous Helmholtz equation with the equivalent current 
density 

j e «(r) = -K^(r)E«(r), (17) 

which justifies the interpretation of Y (r) as an effective 
susceptibility tensor connecting an induced linear polar- 
ization P (r) = £ Y(r)E l (r) of the particle with the 



incident field E 4 (r) . In this formulation one can imme- 
diately recognize that in the case of the resonant scat- 
tering the integral in (|16|l will be defined by the poles of 
the effective susceptibility 



Y(r)= X (r) (Y 



(18) 



as it is required within the linear response theory ap- 
proach (3^ | . 

The integrals in (fTS"]) and (TTB")) are improper at r = r' 
due to a singularity of the scalar Green's function. Con- 
sequently one cannot calculate an effective susceptibil- 
ity without an appropriate regularization. To extract the 
singularity, one can exclude an arbitrary principal vol- 
ume V* around the singular point r, and proceed as in 
(35| to express (fT8|) in the regularized form 



Y (r) = k 2 ( f d 3 r't?(r, y') X 00 e ik ( r '" 

f d 3 r'te(r,r')x(r')e ik ( r '" r )- 
Jv* I 



-(V-[V 3o (r,r')])x(r) 



■^v«x(r)l (19) 



where G(r, r') is the dyadic Green's function 



^(r,r')= ( Y + -V®v) 5 (r,r'), (20) 



and go is the static Green's function 

„A 1 1 



So(r,r') = lim5(r,r 



fc->0 



An r — r' 



(21) 



Yy» denotes the source dyadic 



-in 



av* 



l2 ,R ® n 



(22) 



which accounts for the depolarization of the excluded 
volume V* and depends entirely on the geometry of the 
principal volume, but does not depend on its size or po- 
sition [36|. Here R = r' — r, R — |R|, dV* is a surface 
enclosing principal volume V* and n is its outward unit 
normal. The dyadic Green's function defined in (|20p can 
be calculated explicitly and is equal to 



H(r,r') 



AnR 



ikR 



n Y 



k 2 R 2 

3-3ikR- k 2 R 2 R<S)TL' 



k 2 R 2 



R? 



(23) 
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For points inside the particle , the integral in ([TO)) 
can be regularized in the similar fashion, resulting in 

E(r) = E»+fc 2 ( f dV tf(r, r') V (r') E l (r') 4 

iv* L 

^(V-[V.g (r,r')])^(r)E l (r) 

--^v^(r)E<(r)}. (24) 

For points outside the particle general solution is given 
by 

E(r) = E ,: (r) + k 2 f d 3 r't?(r, r')1? (r') E* (r') . (25) 
Jv 



2.2 Rayleigh approximation 

If the nanoparticle dimensions are much smaller than the 
incident light wavelength, the scattering problem ((SJ) can 
be treated within a quasi-static limit |k| = 0. Letting 
the principal volume to be infinitesimal, the regularized 
integral in (fT9]l can be reduced to [37] 



Xr 

4tt ./,,! 



dV-Jp- = -xX. (26) 



Here the susceptibility of the particle \r is taken out- 
side of the integral, while it is assumed to be constant 
inside the particle. The integration is performed over the 
particle surface dV. In we have introduced the depo- 
larization tensor of the particle L. This tensor depends 
solely on the particle shape. 

Taking into account (I26p one obtains a quasi-static 
effective susceptibility of the small particle in the form 



¥ + X r !> ) 



(27) 



The field induced in the particle is given by the quasi- 
static limit of {2U 



E(r) = E* - t¥ 



E' 



(28) 



which, taking into account that the depolarization tensor 
is symmetric (3(|, can be combined with (|27|) to obtain 
the familiar relation for the electric field in the Rayleigh 
approximation 



E(r) = (? 



X 



E' 



(29) 



Further we calculate the effective quasi-static sus- 
ceptibility for a small metallic nanoparticle. Taking into 
account that the depolarization tensor is symmetric, the 
effective susceptibility (|27p can be written in component 
form as 

= Xr 



1 + XrLi 



(30) 



where 5q are the unit vectors along the main particle 
axes. Assuming for simplicity the Drude model for the 
metal 2 

Xr = , (31) 

with uj p being the plasma frequency and 7 being the 
damping rate, the relation (|30p results in 



•JlLa - uj 2 - i7w 



(32) 



Then the polarization of the metallic particle P = eo^E 1 
for an arbitrary plane harmonic incident field E 4 should 
satisfy the following algebraic equation 

- lj 2 P - i 7 ^P + oj^LuP = eowgifcj <8 x,E l , (33) 

which after taking an inverse Fourier transform results 
in the dynamic equation for a harmonic, damped, driven 
oscillator 

d 2 P dP 



<5<- ' +UJ oP = eoUJpHi^XiE 1 (34) 

with the resonance frequency ujq = lo 2 Lu defined by the 
particle material and shape. Such collective oscillations 
of the electron plasma confined to the metallic particle 
corresponds to a particle plasmon-polariton excitation. 



2.3 Retardation 

With increasing particle size the retardation effects be- 
come important and to obtain an effective susceptibility 
of the particle as well as associated scattered field one 
has to calculate the integrals in fTS"]) and (fTOj) explic- 
itly. To do that, one should choose the principal volume 
carefully. Technically the form and the size of the princi- 
pal volume can be arbitrary [35] . In order to recover the 
quasi-static limit it is convenient to choose the principal 
volume similar in shape to the considered particle. Fur- 
thermore, the principal volume can be chosen infinitesi- 
mally small, resulting in the following exact solution of 
the scattering problem for the points inside the particle: 



E(r) = E i (r) + fc 2 / d 3 
Jv-v* 

tf(r, r') V (r') E* (r') - H^t (r) E* (r) (35) 
with the self-energy tensor given by 

^ (r) = k 2 



3 „' 



dV 



v-v* 



^(r,r')x(r')e ik -( r '- r )-^ X (r). (36) 

Note that in (j3"S")l and ([51?]) L is the depolarization ten- 
sor of the particle and the infinitesimal principal volume 
V* has to have the same shape as the particle under 
consideration. 
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Denoting the result of the integration in (pJB]) by D (r) 
one obtains the effective susceptibility of the particle 



^ (r) = Xr ( 



. L 



-ff(r) 



(37) 



Spatial and frequency dependence of the effective suscep- 
tibility originated from the field retardation is contained 
in the dynamic depolarization tensor (r) . 

To estimate the main effect of the field retardation 
on the effective susceptibility tensor we evaluate the dy- 
namic depolarization tensor. The following approxima- 
tions are used to simplify the analysis. First, we assume 
that the phase difference due to the incident field is neg- 
ligible within the particle, i.e., exp (ik • (V — r)) w 1 in 
(j3"6")h Second, we use only the first few terms in the long- 
wavelength expansion of the dyadic Green's function in 
powers of k to calculate the integral 



f 






g R(g>R 


An 


(- 








A: 2 


(1 


R <g> R 




8tt 




+ R* 



+ 



D7T 



Finally, we use the mean value theorem to approximate 
the resulting integral by its value in the geometrical cen- 
ter of the particle. 

For a centrosymmetric particle, e.g., a sphere, a reg- 
ular cylinder, or a block, the integral of the first term 
in (|38p is equal to zero, the second term results in some 
symmetric shape-dependent tensor k 2 V D ', while the in- 
tegral over the last term results in ik 3 (V/6ir) I . This 
leads to the following form of the dynamic depolariza- 
tion tensor 

V ^ 

()7T 



f)(r) = fcVf)' + ik 3 — I 



(39) 



Taking into account that this dynamic depolarization 
tensor is symmetric, the effective susceptibility (|27[) can 
be written in component form as 

Xr . ^ „ 



l + Xr L tl -k 2 VD' u -ik 3 ^' 



(40) 



which for metallic nanoparticle described by the Drudc 
model (f3~Tj) leads to the following equation for the in- 
duced polarization 



- uj P - i7wP + ujpLii'P + „■ 

K V 



10J 



Gttc 3 



6irc) 
5 x ? -E l 



(41) 



Taking the inverse Fourier transform results in 



d 2 P dP 

_ +7 _ +W0 P_ 



73 



d 3 P 



■ UJ4- 



d 4 P 



75 



d 5 p 

d^~ 



• X,;E l 



(42) 



E 




1 

Figure 2 Definition of the geometrical parameters, radius a 
and length I, of the scattering cylinder as well as the chosen 
body centered coordinate system. Additionally the incident 
angle £ and the polarization basis vector E|| are depicted in 
the incident plane. 
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Figure 3 Comparison of the scattering cross section for 
scattering normal incident light on a gold cylinder (radius 
10 nm, length 100 nm) calculated using the Rayleigh suscep- 
tibility (dashed) and the retarded susceptibility (solid). 



which describes an anharmonic, damped, driven oscil- 
lator. Parameters ujq, W4, 73 and 75 can be easily de- 
duced from the comparison of equations (|41l) and (|42l) . 
As can be clearly seen from equation (|42l) . retardation 
leads (i) to anharmonicity and as a consequence to a 
shift of the resonance frequency (frequency of the par- 
ticle plasmon-polariton) and (ii) to additional radiation 
damping, which broaden the resonance itself. 

In figure [3] the total scattering cross section for nor- 
mal incident light scattered on a regular gold nanocylin- 
der is shown. The geometry of the problem is depicted 
in figure [5J The following parameter are used for calcu- 
lations, particle radius is a = 10 nm and its length is 
I = 100 nm. The relative permittivity of gold is modeled 
by a free electron Drude-Sommerfeld model 



(43) 



with = 9, u^ = 1.36674 • 10 16 s^ 1 and 7 = 7.59297 • 
10 13 s _1 (see |38|). Figure [3] shows both the scattering 
cross section obtained from the effective susceptibility 
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(j2~T)) within the Rayleigh approximation (dashed line) 
and from the retarded effective susceptibility (j4TJ)) . One 
can clearly see that retardation results in a resonance 
frequency shift towards lower frequencies and to a broad- 
ening of the resonance peak. 

2-4 Numerical solutions 

To obtain a more rigorous solution of the scattering 
problem, direct numerical methods have to be used. One 
class of methods discretize the volume integral equation, 
that is equation (|35D with the replacement X (r') E l (r') — > 
E(r'). A popular representative of this class of solvers 
is the discrete dipole approximation (DDA) [39] which 
solves for induced dipole moments of small cubes in 
which the scatterer is decomposed. A numerical method 
also based on an integral solution of Maxwell's equation 
but expressed in terms of surface instead of volume in- 
tegrals is the boundary element method (BEM) jijj. An 
advantage of integral equation based methods compared 
with direct solutions of Maxwell's equations is that only 
the scatterer has to be discretized, in BEM actually just 
its surface, and one does not have to deal with unwanted 
reflections at the borders of the computational domains. 
The approach of solving Maxwell's equation in the fre- 
quency domain directly on a volumetric unstructured 
mesh is called the finite element method (FEM) (4l| . 
The idea in FEM is to approximate the electromagnetic 
field by polynomials of low order on each small mesh vol- 
ume. Plugging this test functions into Maxwell's equa- 
tion together with boundary and connection conditions, 
one yields a system of linear equations to be solved [4lj . 
Similar in both the basic idea and mesh type but formu- 
lated in time domain is the discontinuous Galerkin time 
domain method (DGTD) [42|, Hj|- Another very popular 
time domain method is the finite-difference time-domain 
(FDTD) method Formulated on a structured Carte- 
sian grid and approximating the spatial and time deriva- 
tives by central differences, this algorithm is explicit, i.e., 
no set of linear equations has to be solved. Therefore 
the implementation of the method is relatively simple 
and numerical stability is not a concern as long as space 
and time discretization are connected by the Courant 
relation [44|. On the downside one has to mention that 
one has to choose a very small space discretization to 
accurately describe curved surfaces of scatterers which 
leads to high memory consumption and long computa- 
tion times. 

In figure @] the total scattering cross section calcu- 
lated with different direct numerical methods is depicted 
for scattering of a plane wave under slanting incidence 
(£ = j) on a gold nanowire with radius a = lOnm and 
length I — 200nm (see figure^). Gold is described by the 
equation (I43[) with the associated parameters. The nu- 
merical methods and their implementations are, in par- 
ticular: (i) adda, an open source DDA implementation 
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Figure 4 Scattering cross-section of a gold nanowire (I = 
200 nm, a — 10 nm) under slanting incidence (£ = 21) calcu- 
lated with different rigorous numerical methods. In the top 
panel (a) the first and in the bottom panel (b) the third 
resonance peak are shown. 

[45| | (solid line); (ii) HFSS, a commercial FEM frequency 
domain solver from Ansys Inc.[46j (dashed line); and (iii) 
a self-implemented FDTD code (dotted line). The two 
panels a) and b) show the first and third resonance peak 
respectively. The space discretization in both the DDA 
and FDTD calculation (both methods use a cartesian 
voxel mesh) was set to 1 nm which leads to reasonable 
convergence. The mesh in HFSS is created by an adap- 
tive algorithm so that there is no uniform discretization 
step. As basis functions second order polynomials are 
used. The execution time of DDA and HFSS per fre- 
quency point was about 8 and 6 minutes on one core 
of a modern workstation, respectively. The given execu- 
tion time of HFSS has taken into account the adaptive 
recreation of the computational mesh for each of the 
depicted panels. FDTD, due to its time domain roots, 
yields the whole spectrum in one run which needs about 
250 minutes. As it can be seen in figure 2] the agreement 
between the different numerical method are reasonable 
but not perfect despite high accuracy in sampling and 
long execution times. 

2.5 Semi- analytical methods 

To overcome the dependence on time consuming numer- 
ical simulations we propose some semi-analytical meth- 
ods to calculate the scattering of an incident plane wave 
on a cylindrical wire as depicted in figure [5] Further we 
assume that the radius a of the wire is small compared 
with the wavelength of the incident light (electrically 
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Figure 5 Convergence of the first three resonances (from 
top to bottom panel) of a gold nanowire (I = 200 nm, a = 
10 nm) under slanting incidence (£ = j ) calculated with VC- 
IE method (solid line), SI-IE method (dashed) and Hallen's 
method (dotted). 
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Figure 6 Scattering cross-sections of gold nanowire with 
fixed radius a = 10 nm and different lengths. The full and 
dashed lines show the VC and SI integral equation results, 
respectively. The dotted line is the numerically rigorous ref- 
erence calculated using a DDA method (4t| . 



thin wire). In this case the incident light interacting with 
the wire can be regarded as a function of z alone 



E l (r) a Ege 



II —ikz cos £ 



(44) 



where E denotes the field component along the wire 
axis. First we briefly review how the scattering problem 
is solved in microwave and RF spectral regime. In this 
spectral regime metals can be treated as perfect electric 
conductor (PEC) with the consequence that the induced 
current density has just a longitudinal z-component re- 
siding on the wire interface. Expressed in terms of the 
total current I the induced current density in this case 
can be written as 



j(r) = zl(z 



$(P ~ a) 
2na 



(45) 



By using (|45l) in the Lippman-Schwinger equation ([5]) 
instead of the equivalent current density (U) one derives 
the integro-diflerential equation 



ujjM) 

2tt 



dz 1 



b'g a (ct>',z-z')I(z') (46) 



ii 



where g a (<j)',z — z') = g(r,r') with r = (a, 0,z) and 
r' = (a,cj)',z') in cylindrical coordinates. For a perfect 
electric conductor the tangential electric field at the wire 
interface, E z (a, z), vanishes and (|4T>1) reduces to the well 
known Pocklington integro-diflerential equation (48[. In 
the optical and near-infrared spectral regime metals ex- 
hibit finite conductivity and a skin depth of the order of 
typical dimensions of plasmonic particles. However, as 
long as the wire is thin, regarding the induced current 
as pure surface current should still be a plausible ansatz. 
In contrast, vanishing tangential electric field E z {a,z) 
can not be accepted for wires with finite conductivity. 
So one approach to derive a self-consistent equation to 
determine / is to express the left hand side of (|4"o]) in 
terms of the current / by utilizing the known exact solu- 
tion of the plane wave scattering problem on an infinite 
cylinder [49( . This should work reasonably well for wires 
with high aspect ratio, i.e. I ^> a. The analytical solu- 
tion shows that the electric field inside an infinitely long, 
electrically thin cylinder under plane wave incidence can 
be approximated with high accuracy by 

E(r) a zf(z) Jo (k p p) (47) 

with Jo being the Bessel function of the first kind of 0-th 
order and k p — kyj e r — cos 2 £. From this expression one 
can derive a relation between the current and the surface 
electric field as 



E z (a,z) = Z s I{z) 



(48) 
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with the surface impedance 

„ _ . J (k p a)k p 



2naujeQAe r Ji (k p a) 



(49) 



Relation (|4"5|) can now be used to express the surface 
field on the left-hand side of (|4"r!]) in terms of the current 
/. This yields a self-consistent surface impedance (SI) 
integro-differential equation (Hoi 



Z s I{z) = El nc (z) + i 



dz' 



2ji 



k 2 dz 2 



d<t>'g a {4>',z-z')I(z') (50) 



suitable to calculate the current in thin wires with high 
aspect ratios and finite conductivity. However, a more 
consistent treatment is possible. One can directly use 
(l4"Tj) as an ansatz for the electric field in (|4l))) on both 
sides. The result we obtain is the volume current (VC) 
integro-differential equation 



f(z)J (k p a) = El nc (z) + k 2 Ae r 1 + 



k 2 dz' 



d 3 r' g(a,z;r')f(z')J {k p p') (51) 



which determines the amplitude function f(z). 

Both in equation ([50)1 and in (f5T1) the differential op- 
erator 1 + k~ 2 d 2 can be brou ght inside the integral by 
means of Lee's regularization [35| already used in 
We thus get rid of the necessity to impose additional 
boundary conditions. The regularization step is strictly 
speaking not needed in case of the SI integro-differential 
equation (I50p because it can also be solved in its present 
form by Hallen's approach [48l | using the boundary con- 
ditions I(z) = for z = —1/2 and z = 1/2. However, us- 
ing pure integral equation is the preferred way because 
I(z) is discontinuous at the wire tips so that enforcing 
the current to vanish at the wire ends leads to bad con- 
vergence when solving (|50l) numerically. In contrast the 
VC integro-differential equation can not be solved in its 
present form ([5T1) because we do not know appropriate 
boundary conditions for the electric field at the tips. Af- 
ter regularization both the SI-IE and VC-IE can be dis- 
cretized within a point matching method of moments 
(MoM) scheme [48] . This leads to an n-dimensional ma- 
trix equations in the form 



u = ( r V - U 



E! 



(52) 



where the dimensionality n is the number of slices in 
which the wire is divided. The i-th component of U de- 
notes the unknown function values I{zi) and / (z>) at the 
discrete point Zi for SITE and VC-IE respectively. Like- 
wise the «-th component of E* denotes the z-component 
of the incident field at the slice at position Z{. Obviously, 



r and M have to be calculated differently for SI-IE and 
VC-IE. In both of them, however, r contains the source 
dyadic and all integration over the slice in which the 
singular point resides, whereas the entries of M contain 
the proper integrals over slices without singularity [47| . 
The integrations involved in the calculations of r and 
the entries of M as well as the matrix inversion have to 
be performed numerically. 

To make sure that an appropriate space discretiza- 
tion A (slice thickness) is used when applying the VC-IE 
or SI-IE to realistic scattering problems, their conver- 
gence is studied. Figure [S] shows the dependence of the 
resonance frequencies on the space discretization A for a 
gold nanowire (/ = 200 nm, a = 10 ran) under slanting 
incidence (£ = j). From the top to the bottom panel 
convergence of resonances of increasing order are shown. 
The calculations are performed with (i) the SI-IE solved 
by Hallen's approach (dotted line), (ii) the regularized 
SI-IE (dashed line) and (iii) the VC-IE (solid line). It 
can be seen that (i) the converged result of the SI-IE 
does not depend on the way of solving it and (ii) a dis- 
cretization of A = 1 nm already shows a deviation from 
the converged value of less than ITHz. Thus we choose a 
discretization of A = 1 nm in subsequent studies. In fig- 
ure |B] the total scattering cross section as function of fre- 
quency for gold nanowires of different lengths but fixed 
radius, a — lOnm, under normal incidence (£ = §) are 
shown [47] . The results of VC-IE, SI-IE and numerically 
rigorous DDA method are represented by solid, dashed 
and dotted lines, respectively. One can trace that (i) 
the deviation of both semi-analytical methods from the 
numerical exact DDA reference solution increase with 
decreasing aspect ratios of the wire and (ii) the accu- 
racy of the VC-IE is slightly better than the one of the 
SI-IE. This behaviour was expected because the ansatz 
(|47|) for the internal field taken from the infinitely long 
wire becomes worse for decreasing aspect ratios. It is re- 
markable that for a wire with aspect ratio as low as | 
(I = 50nm) the relative deviation is less than 2%. But 
most important are the differences in execution time. 
The calculation of one frequency point in SI and VC in- 
tegral equation method requires, respectively, about 1 
and 2 seconds on one core of a modern workstation us- 
ing Mathematica [51|. An additional advantage of the 
VC-IE despite its slightly better accuracy in the calcu- 
lation of far field properties like the scattering cross sec- 
tion, is its ability to calculate even the near field of the 
plasmonic wire with high accuracy. To demonstrate this 
a comparison between the internal field obtained with 
the VC-IE and HFSS is shown in figure [7] for a gold 
nanowire (/ = 200nm, a — lOnm) under slanting inci- 
dence (£ = f ). Here the radial component is calculated 
by assuming a solution of the form 



E(p, z) = pE p (p, z) + z/(z) J (k p p) 



(53) 
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Figure 7 Induced electric field at the first (left) and second (right) resonance inside a gold nanowire with radius a — 10 nm 
and length I = 200 nm under slanting incidence (£ = ?, Eg — 1). Each top panel show the z-components on the wire axis and 
the bottom panels the radial p-components 9 nm apart from the axis. The fields are obtained by HFSS (red solid line) and 
VC-IE (black dashed line). 



and enforcing the field to be divergence free. This results 



E P {p,z) 



J\{k p p) d 
k p dz 



/(*)• 



(54) 



The left and right column of figure [7] show the field com- 
ponents for the first and second resonance of the wire 
respectively. The top panel on each column shows the z- 
component of the electric field at the wire axis and the 
bottom panel the p-component in 9nm distance from the 
axis. The p-component on the axis vanishes. The agree- 
ment between VC-IE and HFSS is very good. Thus we 
can conclude that using the calculated internal field in 
the form (|53p , the near field can be calculated with high 
accuracy by means of (|25[) . 



3 Plasmonic dimers 

3.1 Two coupled nanoparticles 

A considerable enhancment of electric field near a reso- 
nant plasmonic nanoparticle opens many application op- 
portunities. Among others, one should mention strong 
interaction between the different type of emission cen- 
ters and the nanoparticle, which results in strong mod- 
ification of emission properties [H^, 53, [H3|. The field 



enhancement also allows to implement strong coupling 
among different nanoparticles themselves. A number of 
intriguing phenomena can be observed in coupled plas- 
monic structures, including Fano-type resonances [55l 
Hfj Ht} , plasmon- induced transparency [HI, [5^ (an ana- 
log of the electromagnetically induced transparency (60| ) 
as well as unusual material properties of metamaterials 

ai. 

To rigorously describe even the simplest type of the 
coupled nanoplasmonic system, i.e., two particles or a 




Figure 8 Sketch of the dimer structure under consideration. 
All relevant model parameters are indicated. 



plasmonic dimer, one has to use direct numerical sim- 
ulations (6ll . IBli . Although a very intuitive coupled os- 
cillator model [62| can explain the basic physical mech- 
anisms of the particle plasmon interaction, it relies on 
phenomenological fitting parameters and is typically ap- 
plicable only to the simplest case of parallel or anti- 
parallel particles. In this section we extend the coupled 
oscillator model to describe the coupling between two 
nanoparticles, which can be arbitrarily oriented with re- 
spect to one another. A sketch of the typical structure 
under consideration is shown in figure [3] 

As was demonstrated in Sec. 12. li the effect of the ex- 
ternal field on a nanoparticle can always be described 
using an effective susceptibility tensor. Introducing par- 
ticle polarizability tensor *cti — eo^t \ (i = 1, 2), the 
polarization (dipole moment) induced in one particle in 
presence of the second one can be written in the form 



( r i) 



< ct, 



E(r0 



(55) 
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with E (rj) being the total electric field at the particle i 
position. Assuming further that electric dipole response 
of the particle can be modeled modeled as a response of 
the point dipole with the same polarizability the total 
field E (rj) can be expressed 



E (r,) = *eti ( E (r<) + — tf (R) £, 



(56) 



with R 



rj, where r^ (rj) is a position of the 



ith (j'th) particle. Here G is the Green's function de- 
termined from equation (I23[) . and the second term is a 
contribution of the equivalent dipole of the second par- 
ticle to the total field at the position of the first particle. 
To find an effective induced polarizability of the parti- 
cles in the dimer, one has to solve the system (|55|) with 
the field expressed by ([551) . For example, substituting an 
expression of the induce dipole moment of the particle j 
in the equation for the dipole moment of the particle i 
we obtain 



Aj - #iE( ri ) 



a 4 H a 4 Kjr ( a. ) a . 



A- 4 



E - 



^W l G > (R)^-^(-R)/i l (57) 



'0 



which can be solved to obtain the unknown effective 
dipole moment fi i 



Mi 



V - ^ G^R)^- G*(-R) 



£0 



r)^, 



Eo — w jEo- 



(58) 



Here the effective polarizability Wj of the particle i in 
the presence of the particle j has been introduced. Please 
note that this expression is valid in the point dipole ap- 
proximation for an arbitrary mutial orientation of the 
nanoparticles. 

A pair of electrical dipoles is general results in effec- 
tive electric dipole, electric quadrupole, and magnetic 
dipole moments. The effective dipole moment /x c ff = 
Ai + A*2 i s given by a simple sum of the individual dipole 
moments /i^ and results in an effective polarizability of 
the plasmonic dimer in the coupled dipole approxima- 
tion given by 



a off 



£ 

i=l,2 



n -1 



a j — G 



(59) 



The electric quadrupole and magnetic dipole moments 
are given respectively by 

j y 

Q = (ri (g + /it! g> ri + r 2 g> /x 2 + A2 ® 1*2) (60) 
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Figure 9 Scattering cross sections of single nanoblocks. 
Lengths are left to right from 60 nm to 120 nm with 10 nm 
step. Width and height of all nanoblocks are 20 nm. 
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Figure 10 Scattering cross section of the gold nanorods 
dimer for left circular polarized light calculated within cou- 
pled dipole model. Solid line is a result of contributions from 
all considered multipole terms. Dashed line is purely electric 
dipole contribution. Dotted line shows combined contribu- 
tion of magnetic dipole and electric quadrupole. Gold rods 
are 100 nm apart and are rotated 45 degree with respect to 
each other. Rods are 80 nm and 75 nm long. 



and 



m = -— (ri x ^ 



r 2 x fi 2 ) 



(61) 



We further use developed couled dipole model to cal- 
culate a scattering cross section of the gold dimer. For 
small particles the dominant contribution to the first 
resonance is dipole in nature. Moreover if a particle is 
elongated, within the spectral range of the longitudi- 
nal plasmon-polariton resonance the contribution of the 
transverse resonances can be neglected. With a good ap- 
proximation corresponding dipole polarizability can be 
described by the Lorentz model with all retardation ef- 
fects included in the renormalized resonance frequency 
and damping (Sec. 12. 3[) . In what follows we model both 
nanoparticles using dipole polarizability of the form 



Ctjflj (g flj, 



(62) 



where the resonance frequency Li j, damping constant 
7j and amplitude fj are determined by the fit the cor- 
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Figure 11 Scattering cross section of the gold nanorods 
dimer for different rod-to-rod distances. Left circular polar- 
ized light is considered. Distance between rods is 100 nm 
(middle panel) and 60 nm (bottom panel). Solid black lines 
are results of the coupled dipole model. Blue lines with dots 
are results of direct numerical simulations performed with fi- 
nite element method (HFSS) [4rj| |. In the top panel scattering 
cross cross sections are shown for individual rods, 80 nm long 
(solid line) and 75 nm long (dashed line). 

responding numerical scattering cross section of the in- 
dividual particles using the Lorentz model. In figure |H] 
scattering cross sections of the gold nanoblocks of differ- 
ent length are shown. Both finite element data (dots) and 
the results obtained using the fitting function (j6"2")) are 
presented, justifying the use of the proposed approach 
to modeling of the individual particles. 

In figure QJJ scattering cross section of the gold dimer 
build from 80 nm and 75 nm long nanoblocks is shown. 
Width and height of both nanoblocks is equal to 20 nm. 
Particles are situated 100 nm apart and at 45° with re- 
spect to each other. Scattering cross section is calculated 
for the left-handed circular polarized light incident per- 
pendicular to the dimer plane, taking into account only 
the electric dipole moment contributions (solid line), the 



v=434 THz v=450 THz 




Figure 12 Far-field directionality diagrams are shown for 
two frequencies corresponding to two maxima in the scat- 
tering cross section. Directionality diagrams are shown in 
the dimer plane. Solid lines are results of the coupled dipole 
model. Lines with dots are results of direct numerical sim- 
ulations performed with finite element method (HFSS) (ir^ . 
Gold rods are 100 nm apart and are rotated 45 degrees with 
respect to each other. Rods are 80 nm and 75 nm long. 



electric quadrupole and magnetic dipole contributions 
(dotted line), and all available multipole contributions 
(dashed line) . One can clearly see that for the consid- 
ered structure the electric dipole contribution is domi- 
nant and contributions from the higher order multipoles 
can be neglected. 

In figures [TT1 and [T2l comparison of the coupled dipole 
model (solid lines) with the direct numerical simulations 
using finite element method (lines with dots) is shown. 
In figure [TJJ scattering cross sections are shown for two 
distances between particles, namely 100 nm and 60 nm. 
One can clearly see that with decreasing the distance 
the quantitative agreement between analytical and nu- 
merical results becomes worse, but the qualitative agree- 
ment remains very reasonable. To appreciate the amount 
of the "cross-pushing" (repulsion) of the original reso- 
nances, the scattering cross sections of the individual 
nanblocks are presented in the top panel of figure [TT] 
Proposed coupled dipole model demonstrate very good 
accuracy at moderate particle separations both for inte- 
gral and differential scattering properties as can be seen 
from figure [T2T where a comparison of the analytical and 
numerical far-field directionality diagrams is shown. 



3.2 Planar chiral meta-atoms 



Having established that coupled dipole approximation 
is valis at least in the case of planar arrangements of 
nanoblocks in the dimer, we proceed to derive the effec- 
tive material parameters. 



Plasmonic nanoparticle monomers and dimers: From nano-antennas to chiral metamaterials 



13 



We begin by rewriting the equation (f5T?f , using . 
in the form 



a e ff = 



1 



acini (H> Mi + ^2/^2 <&> ^2 



1 — OL\aiK 
+ aia 2 n (Ai ® fi2 + A2 ® Ai) 



eff ~ 

=ar/xi 



1 Ai + aj^Aa 



'M2 



a T (Ai ® A2 + A2 ® Ai) 



(63) 



where the unit vectors A? denote the orientation of the 
rods (see figure [5]), andK defines the coupling coefficient 
between the rods: 



so 



Gi (Ai • A2) + Gt 



(Ai • r) ( 



A2 • R 



(64) 



The effective permittivity tensor ^ e ff is derived from 
^c?eff in equation. (|63p as 



e^cff = I — (eoVcoii) 1 ^'eff 



(65) 



In axial representation, ^ e ff can be expressed as 



^cff = T + [af /(2e F cell )] (c+ ® c_ + c_ ® c+) (66) 

where the complex vectors c± = (Ai + ^7± A2 ) determine 
the directions of the optical axes, and rj± are given by 



„eff 



± 



yj(af) 2 -afa% 



voir 



(67) 



The dielectric permittivity tensor ((66)) corresponds 
to an absorbing nonmagnetic crystal, which in general 
has two distinct eigenmodes with different polarizations, 
phase velocities, and absorption coefficients. These eigen- 
modes hi^ and their associated refractive indices 
are given by eigenvectors and eigenvalues of the refrac- 
tive index tensor |63| 



— n e 



x <fe>-l„x 



n 



1/2 



(68) 



Here n is a unit vector of the wave normal. The operator 
n x acts on a vector u in such a way that the result is the 
vector cross product: (n x )u = [n x ul. The tensor A _ 

^ y ^ y ^ y ^ y 

is defined as a pseudoinverse to A (A~ A = I — 
n eg) n) and is introduced in place of the true inverse 

^ y 

tensor ( A _1 ), which does not exist for tensors with zero 
determinant such as n x . 

In figure [T3"h . the real and imaginary parts of the re- 
fractive indices for the two eigenmodes (^1.2) are shown 
for propagation direction orthogonal to the meta-atoms 
plane. One can clearly see that the absorption bands for 
the two eigenmodes are different. This is characteristic 
for dichroic media, and it was shown earlier that circu- 
lar or elliptical dichroism is attributable to the optical 
manifestations of planar chirality [14j|. The eigenvectors 
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Figure 13 (Color online) (a) Refractive indices for two 
co-rotating elliptically polarized eigenwaves hi, 2 of Vh in 
equation (|fc>8[) . (b) Transmission and absorption for incident 
wave with polarization given by hi, 2- (c) Absolute value of 
the complex projection coefficients Ci,2 [see equation (169[) ] 
for RH (solid) and LH (dashed) circularly polarized inci- 
dent wave, (d) Analytically calculated difference in trans- 
mission for LH/RH-polarized incident wave. Here a\ = 13 
mm, <Z2 = 10 mm, d = 10 mm, <f> = 45°, and -0 = 0° 



hi. 2 turn out to correspond to co-rotating elliptical po- 
larizations with orthogonal principal axes (see inset in 
figure I13[) . which is a signature crystallographic prop- 
erty of planar metamaterials (l2l . [l3 |. 

Based on the effective permittivity e e ff) one can de- 
termine the optical spectra using the generalized trans- 
fer matrix approach (see, e.g., [6J|). Figure flUb shows 
the transmission and absorption spectra of the medium 
given by (|66f) for an incident wave with polarization coin- 
cident with the eigenmodes (hi^)- The dips in the trans- 
mission spectra can be seen, in clear correspondence with 
the absorption bands. Due to the small variation in the 
real part of the refractive indices, reflection from such 
an effective medium slab is fairly small and does not 
contribute to the transmission dips. 

For an arbitrarily polarized incident wave, the effec- 
tive medium acts as an absorbing polarizing (dichroic) 
filter splitting the incident field H into two waves with 
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polarizations parallel to the crystal eigenvectors hi. 2 as 

H = Cihi + C 2 h 2 (69) 

where Ci j2 are complex projection coefficients. The trans- 
mission of the effective medium is then determined, on 
the one hand, by the relations between C\^, and on the 
other hand, on the relations between the absorption co- 
efficients for the eigenwaves. 

Figure [T3"b shows the projection coefficients C1.2 for 
circularly polarized incident waves. In this case, it can 
be seen that \C\\ = |C7 2 1. However, in the spectral range 
of strong absorption the coupling of RH and LH inci- 
dent waves to the crystal eigenmodes hi 2 is very dif- 
ferent (C-p 2 7^ Cfa)- Hence, circularly polarized waves 
with different handedness interact with the metamate- 
rial with a different strength, and consequently, have dif- 
ferent transmittance (Tr ^ Tl). For planar geometry, 
reversing the handedness of the incident wave polariza- 
tion is equivalent to reversing the direction of incidence 
or exchange the structure with its enantiomeric counter- 
part (see figure (UJ)b,d). Hence, enantiomeric asymmetry 
in the transmission ATrl = Tr — Tl (figure IT3"H ) can 
be used to quantify the "strength" of planar chiral prop- 
erties in a particular structure. It can be seen that the 
maximum value of ATrl corresponds to the overlap be- 
tween the different absorption bands in the dimer. 

Figure [H] shows the spectra ATrl(oj) for various 
shapes of the dimer, analytical results from the equations 
compared to the results of direct 3D frequency- 
domain numerical simulations |46| . It can be seen that 
the effective medium model offers a good coincidence 
with numerical results. 

When the dimer has an in-plane mirror symmetry, 
the 2D enantiomers become indistinguishable and such 
a structure must necessarily be achiral (ATrl = 0). For 
the geometrical transformation considered, an in-plane 
mirror symmetry is achieved in the following cases: (i) 
for rods of equal length (V-shaped dimer, figure RHb): 
(ii) for parallel rods (Il-shaped dimer, figure [T4"b) : (iii) 
for a T-shaped dimer (figure [JJJi). In all these cases, 
the analytical model correctly predicts the absence of 
enantiomeric asymmetry. 

Otherwise, the dimers are seen to exhibit chiral prop- 
erties, which are stronger when enantiomers are more 
distinct geometrically. So there is an optimum rod mis- 
alignment angle <j> ~ 22.5° in figure [Tib. It is also nec- 
essary that the absorption resonances corresponding to 
individual eigenpolarizations (figure and hence the 
resonances of the individual rods, have some degree of 
spectral overlap. Hence, ATrl depends on the length 
mismatch between the rods in a non-monotonic way: it 
first increases when the dimer deviates from the achiral 
V-shape, reaches a maximum value, and then decreases 
again with a pronounced resonance splitting as Aa be- 
comes greater (see figure [Tib). 

It can also be noticed that ATrl becomes smaller 
as the distance between the rods decreases (figure W£k ) . 



This can be explained by cross-pushing of the absorp- 
tion resonances due to an increase in the inter-rod cou- 
pling for smaller inter-rod distances. The evidence of 
such cross-pushing can be inferred from figure 1111 

If the rods were arranged in a non-planar fashion 
(e.g., as shown in figure [lj:), the magnetic dipole and 
electric quadrupole contribution [equations (1601) and (|61[) ] 
would be expected to contribute to the dimer's electro- 
magnetical response much more significantly. It can be 
shown (the detailed calculation will be available in a 
forthcoming report) that 3D chiral properties and gi- 
ant optical activity effects [10| can be expected from 
such dimers. This is fully supported by the observation 
that some non-planar dimers are geometrically similat 
to known 3D chiral meta-atoms (see figures QJi,c) . 

Note, finally, that we have been considering an iso- 
lated plasmonic meta-atom throughout the paper, while 
real metamaterials contain a multitude of meta-atoms. 
Hence, the results obtained here are directly applicable 
to metamaterials in the approximation that the meta- 
atoms are sufficiently sparse and do not interact with 
each other. Importantly, this also means that chiral prop- 
erties in the metamaterials under study are intrinsic (at- 
tributable to the geometry of the meta-atom itself) rather 
than extrinsic (attributable to the meta-atom arrange- 
ment, as in tilted-cross arrays [(55j|). This agrees with ear- 
lier time-domain simulation results [66|. With the pro- 
posed model extended to include the inter-atom coupling 
(see, e.g., [67]), explicit account of intrinsic vs. extrinsic 
effects in 2D and 3D chiral metamaterials can be given. 
Combining the intrinsic and extrinsic contribution in the 
same metamaterial can be used to maximize its chiral 
properties. 

4 Conclusion 

We have systematically described the basic physics of 
light interaction with plasmonic nanoparticles. Starting 
with the general expression for the polarizability of one 
metallic particle (a plasmonic monomer) in an external 
electromagnetic field, we have explored several ways to 
simplify it for the approximation of the nanoparticle size 
being much smaller than the wavelength (the Rayleigh 
approximation). We have also investigated the correc- 
tions to the Rayleigh approach by accounting for the re- 
tardation effects. We have introduced accurate and nu- 
merically efficient ID semi-analytical method to solve 
the problem of light scattering on metallic nanowires for 
optical frequency range. We have discussed the applica- 
bility of the method and its accuracy in comparison with 
the direct numerical simulations. Applications of these 
results to plasmonic nanoantennas are outlined. 

We have also applied a coupled dipole approach to 
a pair of closely located rod-like nanoparticles (a plas- 
monic dimer). By accounting for the interaction between 
the partices in the dimer, we have arrived at analytic 
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Figure 14 (Color online) Numerical (top) and analytical (bottom) dependence of ATrl on (a) inter-rod distance d, (b) 
difference in rod length a± — <J2, (c) rod misalignment angle <j> as in figure [5] and (c) rod displacement angle ij) for </> — 90° (as 
the dimer changes between T- and L-shaped). Unless specified otherwise, a\ = 13 mm, a,2 = 10 mm, d = 6 mm, <f> = 45°, and 

i> = o°. 



expressions of the dimer's polarizability. The effective 
macroscopic material parameters have then been cal- 
culated. These parameters are attributable to an effec- 
tive medium corresponding to a sparse arrangement of 
nanoparticles, i.e., a metamaterial where plasmonic dimers 
have the function of "meta-atoms". It is shown that pla- 
nar dimers consisting of rod-like particles generally pos- 
sess elliptical dichroism [l4[ and function as atoms for 
planar chiral metamaterials 12, [l3|. It is worth noting 
that we have been able to design "meta-atom" structure 
possessing the previously reported polarization effects, 
while being considerably simpler with respect to analy- 
sis and fabrication. Indeed, it is far easier to fabricate an 
array of nanorods than it is to fabricate an array of more 
complicated objects such as nanogammadions, and the 
difference becomes more pronounced as the length scales 
go down. Hence, the results of this work can form the- 
oretical grounds to new, cheaper, easier-to-make chiral 
metamaterials in the optical domain. 

The authors wish to acknowledge financial support 
from the Deutsche Forschungsgemeinschaft (DFG) Re- 
search Unit FOR 557. 
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